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Abstract. - The utility of a cooperative length scale for describing the dynamics of small 
molecule glass-formers is shown. Molecular Dynamics and Monte Carlo simulations reveal a 
distribution of cooperatively moving fractal events below the temperature Ta at which dynamics 
become caged. Guided by these results, four straightforward methods emerge to recognize Ta 
in experimental data and quantify the length scale that grows on cooling below Ta- This 
length scale is consistent with 4-D NMR experiments which are sensitive to the slow moving 
population. 



Many liquids either cannot crystallize or crystallize sufficiently slowly that they vitrify 
below their glass transition temperature Tg . A fundamental understanding of glass formation 
is still lacking because it has not been firmly established whether the pronounced slowing 
down is simply kinetic in origin or there is an underlying thermodynamic character [1-6]. 
This letter shows that there is a natural length scale for cooperative motion that grows as the 
glass transition is approached. 

The 1965 model of Adam and Gibbs [1] suggests that there should be cooperative motion 
in glass-forming liquids. The size ^ of the cooperative volume is related to the configurational 
entropy of the liquid Sc- At temperatures T, sufficiently above Tg, where the relaxation time 
Ta, and the viscosity r], vary as Tq ^ rj/T ~ exp [E/ksT], all molecules undergo independent 
local Brownian movements without signs of cooperativity. As temperature is lowered, the den- 
sity of the liquid gradually increases and Brownian motion becomes hindered, as neighboring 
particles block each others attempts to move. This crowding leads to cooperative dynam- 
ics [5-9], active for all T below the caging temperature Ta- The onset of cooperativity is also 
accompanied by the observed 'caging effect' in the mean-square displacement of a particle 
between the ballistic and self-diffusive regimes, and a reduction in Sc, causing ^ to grow. 
These changes result in a progressively stronger temperature dependence of rj and Tq at lower 
temperatures. 
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Fig. 1 - Simulation results for the length scale of cooperative motion as a function of reduced 
temperature. From molecular dynamics simulations of Lennard- Jones sphere mixtures (A [11]) and 
from Monte Carlo simulations of the bond fluctuation model at density 0.8 for chains of 10 monomers 
(• [18]). Inset is the temperature dependence of the conflgurational entropy for chains of 10 monomers 
(o [10]). 



Since experimental attempts to identify the length scale for cooperative motion have met 
with Umited success [2-4], the dominant evidence for this quantity is from computer simula- 
tions [8-17]. Simulations have the profound advantages of direct observation of motion and 
straightforward identification of both the size and shape of cooperatively rearranged regions. 
Equilibrium simulations of liquids are not yet possible near Tg, but have been done down to 
0.7Ta [10]. In addition to confirming the essential aspects of the Adam and Gibbs model, sim- 
ulations have provided two novel insights. Instead of a single size scale for cooperative motion, 
there is in fact a broad distribution of size scales below Ta [8, 12-19]. The largest size in this 
distribution f grows rapidly as temperature is lowered, as expected by Adam and Gibbs [1]. 
The second important observation is that the cooperatively rearranging regions are not the 
three-dimensional volumes that were initially proposed, but instead are fractal [8,15-18]. The 
observed fractal dimension (of order 2) clearly shows that the majority of molecules within 
the volume have not participated in the cooperative motion. Consequently, a new model 
for cooperative motion was proposed that accommodates these new insights [20]. All glass- 
forming liquids show a temperature dependence of cooperative size scale, with commensurate 
effects in viscosity and relaxation time, in reasonable accord with the expectation of dynamic 
scaling [20,21] in the temperature range Tc <T < Ta- 

e - T„e-^°/^-^ ^ T^e-^"/'^-^ ^ (T - Tc)-^ (1) 

The critical temperature Tc, is slightly below Tg and is the temperature at which the equilib- 
rium extrapolated values of ^, Tq and rj all diverge. The material-specific activation energy 
apparently depends on whether segmental relaxation (Ea) or viscosity (-E,,) is measured, with 
Ejj/Ea — 1.2 for non-polymeric organic glasses. 

Figure n shows data for ^ as obtained from Molecular Dynamics simulations [8] and new 
results from Monte Carlo (MC) simulations on the bond fluctuation model [10]. In the MC 
simulations the mobile particles are identified as those that move over a time scale of interest. 
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Table I - Temperatures and length scales for five glass-forming liquids, 
al. [28] with the exception PDE where Tg = T(Ta = 100 sec). All Tb are the crossover temperature 
between Vogel-Fulcher fits wtth the exception of D-sorbitol [29] and glycerol [4] from a(3-merging. 

References 
[22,30-32] 
[33-35] 
[30,31,36-41] 
[22,30-32] 
[30,31,40-42] 



Property 


Ta [K] 


Tb [K] 


T, [K] 


Tc [K] 


VvdW [A] 


UTa) [A] 


[A] 


PDE 


355 


325 


294 


278 


4.1 






D-sorbitol 


350 


335 


274 


257 


3.4 


3.8 ±1.8 


27. ± 13. 


OTP 


315 


290 


241 


227 


3.7 


3.5 ±1.0 


52. ± 18. 


salol 


275 


265 


218 


204 


3.5 






glycerol 


280 


262 


190 


173 


2.7 


1.9 ±0.8 


25. ± 11. 



We then look for clusters of these mobile particles and find that their sizes are a function of 
time for short times, but quickly become time independent up to t^. We only consider cluster 
sizes in this intermediate time range. Figurenshows mean sizes of the mobile particle clusters 
as a function of T . Both sets of simulations show that ^ is sensibly independent of temperature 
above Ta, but then grows rapidly when temperature is lowered below Ta- The very different 
nature of the simulations used for the data in Fig. ^ strongly points to the existence of a 
growing length scale below at T < T^. The inset of Fig. ^ shows that 5c as defined by Adam 
and Gibbs also changes character at Ta with a broad crossover between Tb [22] below which 
Eq.^or Vogel-Fulcher should describe dynamics and a much higher temperature above which 
dynamics are Arrhenius [22]. 

Guided by simulations, and owing to the abrupt change in the very nature of relaxation 
at Ta [9, 10,23-25], the caging temperature is easily identified by a variety of experiments 
probing liquid dynamics. We demonstrate this point with four dynamics experiments that 
have broad dynamic range: rotation and translational diffusion of molecular probes, self- 
diffusion, dielectric spectroscopy and rheology. Molecular probe and self-diffusion techniques 
are particulary noteworthy because they provide a model independent measure of Ta- The 
length scales extracted from these experiments, which are in quantitative agreement with 
existing ^(T) data from 4-D NMR [26,27], show a strong temperature dependence only below 
Ta- 

Tc is determined by using literature data on and rj well below Ta for glycerol, o- 
terphenyl (OTP), phenolphthaleinedimethylether (PDE), salol and D-sorbitol and fitting them 
to Eq.n (Tabled). Ta is defined from the crossovers in Fig. El Different dynamic experiments 
provide consistent determinations of both Ta and Tc- Table also reports the crossover 
temperature Tb from Donth [4] where dynamics switch from one Vogel-Fulcher form to another 
[22] . Tb is typically about 20 K below Ta and corresponds to the upper temperature limit 
where dynamic scaling or Vogel-Fulcher quantitatively describe the temperature dependence 
of dynamics. 

The Stokes-Einstein relation expects the rotational relaxation time (t^) and translational 
diffusion coefficient Dt of probe molecules are coupled so that their product is independent 
of temperature. Above Ta the probes diffuse a distance of order their own size in the time it 
takes for the probe to rotate. However, at Ta these two dynamics decouple, a fact which can 
be used to establish a length scale. To understand this, consider a bimodal distribution of (j)f 
fast and (j)s — I — (f)f slow particles. The average rotational time is dominated by the slow 
particles 

(Tr) = (f>sTs + (f>fTf = (l)sTs (2) 

since ips ^ (j^f ^-iid Ts ^ Tf- The diffusion coefficient is the sum of fast and slow contributions 
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Fig. 2 - Dynamic data for OTP, PDE and salol above and below Ta- Temperature dependence of 
viscosity (□), Ta (o), probe diffusion (x) and self-diffusion (+). Solid curves are fits of Eq.0to data 
below Ta and Arrhenius fits to data above Ta for both viscosity and dielectric data. The dashed line 
is Ta taken from the extrapolation of dynamic scaling to the Arrhenius temperature dependence of 
viscosity from high-T. See Table Q for references. 



Dt = cf>sDs + cl^fDf = ^ses/QTs + h^f/Q^f (3) 

where = SDs^s and = 6DfTf. Since (j)s ~ 1, the product 6Dt{Tr) ^ ^^^^ + (pjCps^'^fTs/Tj 
can be used to define a length scale, 

i^^s^ VmM. (4) 

Fig. OlA shows the temperature dependence of this length scale for OTP and it quantitatively 
agrees with 4-D NMR, which is known to target the slow contribution Hence, ^ 
(j)f<j)s^fTs/Tf, justifying Eq.H 

In the calculation of ^, (r^) was interpolated [43] using the temperature dependence of 
viscosity from various sources [37,40,41]. In OTP, ^ of anthracene (filled circles in Fig.|2lA) 
shows a strong temperature dependence over the entire experimental temperature range be- 
cause all T < Ta, whereas the larger tetracene probe (filled triangles in Fig. El A) is clearly 
temperature dependent below Ta and independent of temperature above Ta- 




Fig. 3 - The temperature dependence of the cooperative length of OTP, glycerol and D-sorbitol. 
Filled symbols denote absolute measures of ^ from 4-D NMR (■ with error bars [27]) and calculations 
using probe dynamics (• anthracene [36]; ▲ tetracene [36]). Open symbols denote measures of the 
temperature dependence of ^ from self-diffusion (o [30,38]) and dielectric measurements (A [31,34]) 
which have been vertically shifted into agreement with the absolute measures of ^. Dashed lines 
denote Ta and Tg (values in Table|IJ. The solid line is the slope of —3/2 expected by dynamic scaling 
(Eq.Cl. 



Ta is the temperature below which ^ becomes temperature dependent. Cooperativity 
causes the breakdown of the Stokes-Einstein relation, as is made evident when Dtrj/T becomes 
temperature dependent at temperatures below Ta (Fig. EI- 

Ideally a probe molecule would have the same shape, size, polarity and properties as the 
liquid matrix in which it was inserted. ^H-NMR can provide measurements of translational and 
rotational se/f-diffusion coefficients [38] . This allows for direct measurement of ^ using Eq. 01 
Although measurements of (r^) exist [38], they are rare. In the case where measurements 
of (Tr) are lacking, we assume that rj/T properly describes the temperature dependence of 
(r^) [38]. This assumption allows for measurements of Ta that stand in agreement with the 
other measurements of T^, but with the limitation of not being able to give an absolute 
measure of ^. Each of the glass-formers in Figure |21 show that Ta measured using dielectric 
and viscosity data occurs at the same temperature as the breakdown in the Stokes-Einstein 
relation for self-diffusion standing in agreement with the small probe data. 
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In Fig. ISlA, the ^ of OTP from self-diffusion measurements are plotted along with the 
absolute measurements of ^ provided by 4-D NMR {^{252K) = 2.9 ± 1 nm, [26]). While each 
technique measures ^ differently, they all stand in quantitative agreement. The ^ from self- 
diffusion data also shows qualitative agreement with the molecular probe data. While Ta has 
been determined by applying dynamic scaling to viscosity and dielectric measurements, results 
from probe molecules and self-diffusion data are in agreement since ^ becomes temperature- 
independent above Ta- 

Equation ^ can be used to estimate ^ from 77 or r^, since Tc, and Ejj have been 
previously established by fitting Eq. ^ to data. 

^ ~ ^l/6g-Sc/6fcBT Tc<T <Ta (5) 

In Figs. OlB and ElC, dielectric data have been used to plot ^ of glycerol and D-sorbitol 
along with measurements from 4-D NMR and in the case of glycerol, self-diffusion. ^ from 
dielectric and self-diffusion data were vertically shifted into agreement with ^ from the absolute 
measurements that have been provided by 4-D NMR. The resulting data obey the slope of 
—3/2 expected by dynamic scaling (Eq. The ^ plotted in Figures |3lB and|3|C shows a 
clear temperature dependence below Ta- 

£,{Ta) and S,{Tg) are determined from Fig. O and listed in Table HJ Uncertainty in the 
measurement of ^ by 4-D NMR is extrapolated to Ta and Tg in cases where 4-D NMR provides 
the only absolute measure of ^ used in the determination of ^, as is the case with glycerol and 
D-sorbitol. The van der Waals sphere radii r^dw of each of the glasses was calculated from 
atomic radii using the procedures of Edward [44] and are included in Table HI The Tvdw for 
D-sorbitol, OTP and glycerol are all within the calculated range of ^(Ta), suggesting that the 
magnitudes of cooperative size calculated herein are reasonable. 

Four robust experimental methods for determining the caging temperature have been 
identified. Of these methods, estimation of the length scale for cooperative motion is best 
done from measurements of probe diffusion and rotation. However, this has only been done 
for a select few glass-formers. Far more convenient techniques of dielectric spectroscopy and 
rheology can determine the temperature dependence of the length scale. 

Despite the fact that the cooperative volume is fractal instead of space-filling, the essential 
features of the Adam-Gibbs model [1] are correct. The cooperative size does indeed grow 
rapidly as temperature is lowered below the caging temperature. Over a temperature range 
extending from Ta to Tg, ^ has been observed to increase by an order of magnitude from the 
van der Waals radius of each molecule, making the insight of Adam and Gibbs particularly 
noteworthy. 

We thank the National Science Foundation (DMR-9977928 and DMR-0422079) for funding. 
REFERENCES 

[1] Adam G. and Gibbs J., J. Chem- Phys, 43 (1965) 139. 
[2] SiLLESCU H., J. Non-Cryst. Solids, 243 (1999) 81. 
[3] Ediger M., Annu. Rev- Phys. Chem,-, 51 (2000) 99. 

[4] DONTH E., The glass transition: relaxation dynamics in liquids and disordered materials 

(Springer, New York) 2001. 
[5] KOB W., J. Phys. Condens. Matt, 11 (1999) R85. 

[6] Angell C, Ngai K., McKenna G., McMillan P. and Martin S., J. Appl. Phys., 88 (2000) 
3113. 

[7] KiSLiUK A., Mathers R. and Sokolov A., J. Polym. Set. B: Poly. Phys., 38 (2000) 2785. 
[8] Glotzer S., J. Non-Cryst. Solids, 274 (2000) 342. 



B. M. Erwin et ai: ENHANCED COOPERATIVITY IN GLASS-FORMING LIQUIDS 



7 



[9] Binder K., Baschnagel J. and Wolfgang P., Prog. Polym. Sci., 28 (2003) 115. 
[10] Kamath S., Colby R. and Kumar S., J. Chem. Phys., 116 (2002) 865. 
[11] Lacevic N., Starr F., Schr0der T., Novikov V. and Glotzer S., Phys. Rev. E, 66 (2002) 
030101. 

[12] Johnson G., Melcuk A., Gould H., Klein W. and Mountain R., Phys. Rev. E, 57 (1998) 
5707. 

[13] Glotzer S. and Donati C., J. Phys: Condens. Matter, 11 (1999) A285. 
[14] GiOVAMBATTiSTA N., BuLDYREV S., Starr F. and STANLEY H., Phys. Rev. Lett., 90 (2003) 
085506. 

[15] Donati C., Douglas J., Kob W., Plimpton S., Poole P. and Glotzer S., Phys. Rev. 

Lett, 80 (1998) 2338. 

[16] Donati C., Glotzer S., Poole P., Kob W. and Plimpton S., Phys. Rev. E, 60 (1999) 3107. 
[17] Muranaka T., Prog. Theor. Phys. Supp., 138 (2000) 217. 
[18] Kamath S., Doctoral dissertation (Pennsylvania State University) 2003. 
[19] Mel'cuk a., Ramos R., Gould H., Klein W. and Mountain R., Phys. Rev. Lett, 75 (1995) 
2522. 

[20] Colby R., Phys. Rev. E, 61 (2000) 1783. 

[21] Erwin B. and Colby R., J. Non-Cryst Solids, 307-310 (2002) 225. 
[22] Stickel F., Fischer E. and Richert R., J. Chem. Phys., 105 (1996) 2043. 
[23] Wolfgardt M., Baschnagel J., Paul W. and Binder K., Phys. Rev. E, 54 (1996) 1535. 
[24] Baschnagel J., Wolfgardt M., Paul W. and Binder K., J. Res. Nat. Inst. Stand. Tech., 
102 (1997) 159. 

[25] Zhang W., Zou X., Jin A. and Tian D., Phys. Rev. E, 61 (2000) 2805. 

[26] Reinsberg S., Heuer A., Doliwa B., Zimmermann H. and Spiess H., J. Non-Cryst. Solids, 

307-310 (2002) 208. 
[27] Qiu X. and Ediger M., J. Phys. Chem. B, 107 (2002) 459. 

[28] BOHMER R., Ngai K., Angell C. and Plazek D., J. Chem. Phys., 99 (1993) 4201. 

[29] Wagner H. and Riciiert R., J. Non-Cryst Solids, 242 (1998) 19. 
[30] Chang I. and Sillescu H., J. Phys. Chem. B, 101 (1997) 8794. 
[31] Stickel F., Doctoral dissertation (Mainz University) 1995. 
[32] Heuberger G. and Sillescu H., J. Chem. Phys., 100 (1996) 15255. 

[33] Nakheli a., Eljazouli A., Elmorabit M., Ballouki E., Fornazero J. and Huck J., J. 

Phys: Condens. Matter, 11 (1999) 7977. 
[34] NozAKi R., Suzuki D., Ozawa S. and Shiozaki Y., J. Non-Cryst Solids, 235-237 (1998) 
393. 

[35] Naoki M. and Kashima S., J. Phys. Chem., 97 (1993) 12356. 

[36] Cicerone M., Blackburn F. and Ediger M., J. Chem. Phys., 102 (1995) 471. 

[37] Plazek D., Bero C. and Chay I., J. Non-Cryst Solids, 172-174 (1994) 181. 

[38] FujARA F., Geil B., Sillescu H. and Fleischer G., Z. Phys. B, 88 (1992) 195. 

[39] Greet R. and Turnbull D., J. Chem. Phys., 46 (1967) 1243. 

[40] CUKIERMAN M. and Uhlmann D., J. Non-Cryst Solids, 12 (1973) 199. 

[41] Laughlin W. and Uhlmann D., J. Phys. Chem., 76 (1972) 2317. 

[42] SCHROTER K. and Donth E., J. Chem. Phys., 113 (2000) 9101. 

[43] Bainbridge D. and Ediger M., Rheol. Acta, 36 (1997) 209. 

[44] Edward J., J. Chem. Educ, 47 (1970) 261. 



